Testing Small CPAS Parachutes Using HIVAS 
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The High Velocity Airflow System (HIVAS) facility at the Naval Air Warfare Center 
(NAWC) at China Lake was successfully used as an alternative to flight test to determine 
parachute drag performance of two small Capsule Parachute Assembly System (CPAS) 
canopies. A similar parachute with known performance was also tested as a control. Real- 
time computations of drag coefficient were unrealistically low. This is because HIVAS 
produces a non-uniform flow which rapidly decays from a high central core flow. Additional 
calibration runs were performed to characterize this flow assuming radial symmetry from 
the centerline. The flow field was used to post-process effective flow velocities at each 
throttle setting and parachute diameter using the definition of the momentum flux factor. 
Because one parachute had significant oscillations, additional calculations were required to 
estimate the projected flow at off-axis angles. The resulting drag data from HIVAS 
compared favorably to previously estimated parachute performance based on scaled data 
from analogous CPAS parachutes. The data will improve drag area distributions in the next 
version of the CPAS Model Memo. 
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u = Local flow velocity 

V air = Total airspeed relative to air mass 

I. Introduction 

O NE of the challenges in the development of the Capsule Parachute Assembly System (CPAS) is testing an 
evolving design with a limited budget. Building on two previous generations of testing, CPAS is currently 
conducting Engineering Design Unit (EDU) tests with the full Orion architecture as shown in Figure 1. The EDU 
flight test series has successfully obtained data from instrumented Drogue and Main parachutes. However, the 
available test techniques 1 are limited in altitude, making it difficult to test the Forward Bay Cover Parachutes 
(FBCPs) in flight-like conditions. While the mortar-deployed Pilot parachutes have been demonstrated during 
Cluster Development Tests (CDTs), instrumentation to obtain quantitative data has only achieved limited results. 
Further, CPAS cannot use data obtained from first generation (Gen I) Pilot Development Test (PDT) series 2 * because 
the design changed from a ringslot to a conical ribbon. Pilot parachutes were not used during Gen II testing. 
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Figure 1. Orion parachute deployment sequence. This paper will focus on the highlighted parachutes. 


The High Velocity Airflow System (HIVAS) at the Naval Air Warfare Center (NAWC) at China Lake, 
California provides a less expensive alternative to test steady-state parachute performance than flight testing. The 
concept is to create a high-speed flow with a series 
of TF-33 PI 1 engines ducted together. There are two 
HIVAS assemblies at NAWC China Lake. This test 
used the assembly with four engines, as shown in 
Figure 2. This same unit was used to test a candidate 
drogue system for the X-37 approach and landing 
test vehicle (ALTV). 3,4 ’ 5 HIVAS allows using the 
full scale hardware while a wind tunnel would 
require either a large cross section (and therefore 
expensive run time) or using subscale models with 
inherent differences in stiffness and resolution. 

However, because the deployment is not 
representative of flight, inflation parameters cannot 
be determined from HIVAS. 
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Because of the low HIVAS test cost, CPAS decided to utilize the facility to test both FBCP and Pilot parachutes. 
In addition, an Air Force Subscale Aerial Target (AFSAT) Drogue parachute was tested as a control because it has 
known performance. 6 All three parachutes are manufactured by Airborne Systems North America. 


II. Test Procedure 

The test was intended to measure parachute performance at a series of desired dynamic pressures. To measure 
the development of flow velocity, V air , a flow sensor is located just past the HIVAS nozzle exit, and a flow sensor is 
located at the approximate location where the specimen parachute will be when inflated as shown in the top section 
of Figure 3. These velocities are converted to dynamic pressure according to Eq. 1, where p is the ambient density. 

q_ 2' p ' Vair (i) 

The flow sensor at the specimen location must be removed when the parachute is installed as shown in the 
bottom section of Figure 3. In order to achieve the desired flow at the parachute, the nozzle flow reading is mapped 
to specimen flow readings via calibration runs. During each calibration run, the HIVAS throttle is adjusted until the 
specimen sensor reads each desired dynamic pressure. At each setting, the corresponding nozzle dynamic pressure 
reading is noted. When the test parachute is installed, the previously-determined nozzle readings are dialed-in to 
provide the desired flow at the specimen. Calibration runs were performed before and between each parachute 
installation, to account for deviations from the nominal. Ambient atmospheric readings are also recorded, to 
compensate for the rapid heating of the desert environment. 


Calibration Configuration 




Figure 3. The flow is calibrated with a flow sensor (top). The block with the sensor is removed when the 

parachute is attached to a strain link (bottom). 


Velocity and dynamic pressure histories of the nozzle and specimen from the first calibration run are shown in 
Figure 4. 
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Desired qbar values at 
specimen: 



Figure 4. Measured velocity (left) and dynamic pressure (right) at nozzle and specimen. 


Because the flow is non-uniform, additional calibration runs were performed at locations approximately along 
the skirt of each parachute, shown in Figure 5. The ratio of projected diameter to reference diameter, D p /D 0 , was 
assumed to be 0.70 based on Table 5-2 of Ref. 7. Later analysis of photographs confirmed this result. Two of the 
flow surveys were inadvertently performed about one foot below the nozzle centerline, but this oversight was 
minimized by assuming flow is proportional to the radial distance from the centerline for analysis. 
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Figure 5. The flow survey locations were intended to be along the centerline and skirt edges. The initial 
survey was inadvertently about 1 ft below the centerline. 


The assumption of radial symmetry allows for polynomials to be fit to the four calibration points. A second order 
polynomial curve fit was used to characterize the local flow, u, as a function of radius, r. The dynamic pressure and 
velocity profiles are plotted in Figure 6. 
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Figure 6. Dynamic pressure and velocity settings at the nozzle exit (left) and corresponding profiles at 

parachute location (right). 


III. Momentum Flux Correction Factor 

Parachute performance is measured as drag coefficient where the measured force, F, is normalized by the 
dynamic pressure and reference area, S 0 , as in Eq. 2. 
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Figure 7. Non-uniform flow and equivalent flow velocity. 


However, because the measured flow 
velocity is non-uniform, the V air term must be 
computed over the given parachute inlet area. 
The equivalent V air will have the same mass 
flow rate, but a lower momentum flux than the 
actual non-uniform flow. A momentum flux 
correction factor, /?, allows the use of an 
equivalent velocity corresponding to the given 
mass flow rate. The factor is defined as area- 
averaged integral of the square of the 
normalized local flow velocity, according to 
Eq. 3. 8 Therefore, the correction factor needs to 
be computed for the test data at each throttle 
setting for each parachute. Two different 
integration methods were used: a polar 

coordinate formulation and a Cartesian 
coordinate formulation. 



( 3 ) 


A. Polar Coordinate Formulation 

For cases where the parachute is centered along the flow axis, 
polar coordinates are the most convenient. Figure 8 shows the 
polar coordinates used to evaluate the integral. The differential area 
is defined from a differential radius as in Eq. 4. Substituting the 
differential area and the area of a circle into Eq. 3 results in Eq. 5. 
The expression simplifies to Eq. 6. Although this integral could be 
solved analytically by substituting each flow field polynomial 
curve fit, it was simpler to compute each factor with numerical 
integration. 


d A « 2 • n • r • dr 


( 4 ) 
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n 



( 5 ) 



Figure 8. Differential flow area using polar 
coordinates. 
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The resulting equivalent velocity for both size parachutes is compared with the flow field at the penultimate 
throttle setting in Figure 9. Because the FBCP is the smallest parachute, its equivalent airspeed is higher than that for 
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the larger parachutes. This is because the smaller parachute sees more of a concentrated central core flow. However, 
this result is only valid during the brief times when the FBCP is aligned with the centerline. 
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Figure 9. Resulting non-uniform flow and equivalent flow velocity using polar coordinates. 


B. Cartesian Coordinate Formulation 

The polar form of this equation is only useful if the 
parachute is centered along the nozzle centerline. However, the 
FBCP was noticeably less stable than the other parachutes, and 

oscillated about the centerline. To account for a given off-axis J 

angle, 0 \ the momentum flux was also computed using 

Cartesian coordinates. The off-axis angle is shown in Figure 

10. In this case, the differential area is the differential width 

multiplied by differential height, as in Eq. 7, and / 3 is computed 

according to Eq. 8 

dA « dy • dz (7) 



To evaluate this integral, a Cartesian mesh is established and the equation is applied only at points where the 
skirt is projected. This is illustrated in Figure 11. 
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Figure 11. Sample Cartesian meshes to evaluate momentum flux correction factor. 


Sample projected flows are shown for four offset angles in Figure 12. The equivalent flows at the (9=0° 
condition were similar to the results using the polar coordinate system. As the distance moves away from the 
centerline, the velocity curve fits will eventually have negative flow rates, which are not physically possible. 
Therefore, any reversed flows are bounded at zero. 
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Figure 12. Resulting non-uniform projected flows at various off-axis angles using Cartesian coordinates. 
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IV. Load Measurements 


Parachute load measurements were taken by a strain link rated to 10k lbf with approximately 10 lbf resolution. 
Data were recorded at 1 00 Hz. 


A. AFSAT Drogue Parachute 

The AFSAT Drogue parachute was the first parachute tested because it had known performance. The flow at the 
nozzle and the resulting loads trace is shown in Figure 13. 

Parachute loads 

Nozzle qbar settings corresponding to desired 

corresponding to desired S P ec i men c lbar values 



Figure 13. Measured nozzle dynamic pressure (left) and AFSAT Drogue load (right). 


B. EDU Forward Bay Cover Parachute 

Post-test inspection revealed that the confluence keeper on the FBCP failed during the test. This is not a 
significant concern for the operational use of the FBCP in flight because the length of time during which the FBCP 
will be active in flight is much shorter than the duration of the test and the failure mode is not detrimental to the 
purpose of the parachute. Increasing the suspension line convergence angle will tend to increase the projected area 
and drag force. 



Figure 14. FBCP load trace. 
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The FBCP oscillated significantly in all 
radial directions more than the other 
parachutes during the test. Several methods 
were considered to account for the 
oscillation in the drag coefficient 
calculation. The JSC-KX Image Science & 

Analysis Group tracked the FBCP vent 
center coordinates in a longitudinal video. 

A small amount of tracking data is 
superimposed in Figure 15, demonstrating 
the FBCP moved back-and-forth with some 
coning motion. 

Correlating the load trace with the . 

instantaneous off-angle observed in video is ; 

problematic due to the slow frame rate (30 
frames per second) and low resolution of 
the standard definition (SD) video. Also, 
while the video is time-coded, the load data 
did not have an absolute time stamp, 
making synchronization difficult. 

A dual-exposure image in Figure 16 shows the extent of the oscillation, estimated to be about five degrees from 
photographic analysis. Because the flow velocity tapers off from the centerline, the effective flow velocity at this 
position is significantly less than the core flow, as shown in Figure 17. 
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Figure 17. FBCP effective flow velocity comparison of axial flow and largest off-axis angle. 


Another method considered was to 
assume that the maximum loads would be 
encountered when the FBCP was along the 
centerline and the minimum loads would 
occur at the farthest off-axis angle. However, 

a close inspection of the video and data - 

indicate that the observed oscillation * 

frequency was considerably faster than the 

variation in the load trace. A typical J 

oscillation period of about one second is 

shown in Figure 18. However, the tracking 

data showed that the FBCP had oscillated 

about ten times in that period. Only small 

load variations on the order of the 

measurement uncertainty occurred during a 

typical visible cycle. Therefore, variations in 

the loads data are most likely caused by 7 8 

perturbations in the HIV AS throttle setting 
and flow output. 

Knowing the data varied with time, it was 

desired to compute a flow corresponding to the “average” off-axis angle without the efforts and uncertainty of a 
time-varying video analysis. The FBCP oscillation can be modeled as a repeating sinusoidal function. The absolute 
amplitude was assumed to be five degrees, based on the side-view photograph. The time average value for half- 
period of a sinusoid function is the amplitude multiplied by 2/n. Therefore, the average flow velocity was computed 
for an off-axis angle is 5°-2 In « 3.18°. Assuming a large number of periods occur during the averaging window, this 
result is independent of the oscillation frequency. The effective flow at this angle was then used with the time- 
averaged load data to compute drag coefficient. 

C. EDU Pilot Parachute 

As with the FBCP, the confluence keeper on the Pilot parachute failed during the test. Again, this was not a 
significant concern to the program for the same reasons. The load trace is shown in Figure 19. The confluence 
failure is shown on video captures and noted on the loads trace. 
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Confluence 



Figure 19. Pilot parachute load trace. 


V. Drag Coefficient Calculations. 

Real-time drag coefficient calculations during the test did not consider the non-uniform flow and the results were 
unreasonably low. Post-test analysis was necessary to properly normalize the drag with the momentum-flux 
corrected velocity. 

A further correction was applied to the load data to account for the manual adjustment in the HIVAS throttle 
setting. A scaling factor was computed at each throttle setting by dividing the measured nozzle dynamic pressure 
with the corresponding nozzle dynamic pressure during the calibration runs. This assumes a linear relationship 
between the dynamic pressure at the nozzle and the dynamic pressure at the specimen. 

A. AFSAT Drogue Parachute 

The AFSAT results are summarized in Table 1. Each of the four throttle settings are listed in a separate column. 
The readings before flow corrections are shown in red. Because the canopy was relatively stable, the airspeed 
corrections were made assuming axial flow using the polar coordinate formulation. Scaling the loads based on the 
measured nozzle readings resulted in only a small adjustment to each mean load. 

The computed drag coefficients are notably higher at lower throttle settings than the later results. Airborne 
Systems provided a drag coefficient of about 0.55 from previous AFSAT Drogue test data. Therefore, the computed 
Drogue coefficient from HIVAS of 0.606 is about 10% higher than expected. 


Table 1. AFSAT Drogue Parachute Results with Axial Flow 


Mean Airspeed, V a i r (knot) 

131.3 

165.3 

214.8 

250.1 

Corrected Airspeed, V air (knot) 

93.0 

121.1 

172.5 

200.8 

Corrected Airspeed, V air (ft/s) 

157.0 

204.4 

291.1 

338.9 

Mean qbar (lb/ft 2 ) 

51.9 

81.7 

135.3 

180.0 

Corrected qbar (lb/ft 2 ) 

25.6 

43.1 

86.2 

115.2 

Mean Load, F (lbf) 

1,295 

2,088 

3,578 

4,836 

Scaled Mean Load, F (lbf) 

1,268 

2,139 

3,701 

4,901 

Drag Coefficient, C D 

0.649 

0.652 

0.563 

0.558 

Average Drag Coefficient, C D 

0.606 


B. EDU Forward Bay Cover Parachute 

The initial computation of the FBCP drag coefficient assumed the axial flow formulation, with results shown in 
Table 2. The last two throttle settings occurred after the confluence failure. One would expect that the loss of the 
confluence would slightly open the canopy, resulting in a higher drag coefficient, but the computed drag coefficient 
actually dropped. 
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Table 2. EDU Forward Bay Cover Parachute Results with Axial Flow 


Mean Airspeed, V a i r (knot) 

131.3 

165.3 

214.8 

250.1 

Corrected Airspeed, V a i r (knot) 

121.8 

151.2 

207.1 

232.4 

Corrected Airspeed, V air (ft/s) 

205.5 

255.1 

349.6 

392.2 

Mean qbar (lb/ft 2 ) 

51.9 

81.7 

135.3 

180.0 

Corrected qbar (lb/ft 2 ) 

43.7 

66.7 

123.2 

152.8 

Mean Load, F (lbf) 

889 

1,409 

2,309 

3,072 

Scaled Mean Load, F (lbf) 

895 

1,519 

2,413 

3,076 

Drag Coefficient, C D 

0.532 

0.592 

0.509* 

0.523* 

Average Drag Coefficient, C D 

0.539 


*After confluence failure 


The drag coefficient was also computed using the previously mentioned average off-axis angle of 3.18° with 
results listed in Table 3. Because the effective flow is reduced, the drag coefficient is increased for the same 
measured drag force. 


Table 3. EDU Forward Bay Cover Parachute Results with Mean Off-Axis Flow 


Corrected Airspeed, V air (knot) 

110.0 

138.5 

192.0 

219.2 

Corrected Airspeed, V air (ft/s) 

185.7 

233.8 

324.1 

369.9 

Corrected qbar (lb/ft 2 ) 

35.7 

56.0 

105.9 

135.9 

Scaled Mean Load, F (lbf) 

895 

1,519 

2,413 

3,076 

Drag Coefficient, C D 

0.652 

0.705 

0.592* 

0.588* 

Average Drag Coefficient, C D 

0.634 


* After confluence failure 


The drag performance of all CPAS parachutes are documented as drag area probability distributions based on 
flight test data in the “Engineering Development Unit Operating Modeling Parameters Version 11” (aka the “CPAS 
Model Memo vll”). 9 These distributions are explained in Ref. 10. However, the FBCPs had not yet been flight 
tested as of the current model release. The FBCP design is similar to the CPAS Drogue design and both types of 
parachutes have a suspension line length ratio (L s /D 0 ) of 2.0. Therefore, FBCP drag area distribution was determined 
by scaling full open EDU Drogue drag area data according to Eq. 9. The scaled Drogue data were used to define the 
drag area probability distribution published in Model Memo vll. 


<r q\ -tr q\ (Sq)fbcp _ ™ (Dq)fbcp 

VMZfvFBCP v'^D^/Drogue’ x vMD^vDrogue* . .2 

Wo /Drogue v^o /Drogue 


( 9 ) 


The distribution of scaled EDU Drogue test data is shown in blue in Figure 20. The HIVAS data from both 
methods have been appended to the histogram to evaluate how including HIVAS data would affect the distribution. 
The axial method FBCP results were generally lower than the scaled Drogue data and the off-axis results were 
higher. Due to uncertainties in the test data, a decision was made to include data from all three sources to generate 
the best fit normal (Gaussian) distribution. This results in a distribution most likely to encompass the true 
performance. The “nominal” value is the median of all the data. Because the FBCP has only half the number of 
gores as a CPAS Drogue (12 vs. 24), the FBCP was expected to be less efficient than the scaled Drogue data. 
However, the scaled Drogue data were obtained in the presence of a test vehicle forebody drag, while HIVAS flow 
is not significantly obstructed. These effects may have counteracted each other such that the composite distribution 
has a similar mean value to the published distribution. The addition of the HIVAS data has widened the distribution, 
increasing the standard deviation by about 35%. This distribution will be revised as planned flight test data becomes 
available. 
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Figure 20. Comparison of expected FBCP drag area with HIVAS results. 


C. EDU Pilot Parachute 

The Pilot parachute did not exhibit significant oscillations, so the drag coefficients listed in Table 4 were 
computed using the axial flow method. The drag coefficient fell with increasing flow and then rose after the 
confluence failed (the last column). 


Table 4. EDU Pilot Parachute Results with Axial Flow 


Mean Airspeed, V a i r (knot) 

131.3 

165.3 

214.8 

250.1 

Corrected Airspeed, V air (knot) 

93.0 

121.1 

172.5 

200.8 

Corrected Airspeed, V a i r (ft/s) 

157.0 

204.4 

291.1 

338.9 

Mean qbar (lb/ft 2 ) 

51.9 

81.7 

135.3 

180.0 

Corrected qbar (lb/ft 2 ) 

25.3 

42.8 

85.4 

113.8 

Mean Load, F (lbf) 

1,326 

1,919 

3,176 

4,939 

Scaled Mean Load, F (lbf) 

1,317 

1,816 

2,999 

4,835 

Drag Coefficient, C D 

0.687 

0.589 

0.488 

0.569* 

Average Drag Coefficient, C D 

0.583 


* After confluence failure 


The EDU Pilot parachutes have been successfully used in all EDU flight tests. However, only limited data has 
been acquired and direct measurements of drag area may not be possible in flight because the each Pilot parachute 
becomes unloaded as its corresponding Main parachute deploys. The Pilot parachute conical ribbon design is similar 
to a smaller version of the CPAS Drogue design. The EDU Pilot design has an L s /D 0 ratio of 1.15. Therefore, the 
Pilot drag area distribution was determined by scaling full open CPAS Gen I & II Drogue data (L s /D 0 =1.5) rather 
than the EDU Drogue data (L s /D 0 = 2.0). The data were scaled according to Eq. 10. 






( 10 ) 
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Only a limited set of four clean- wake Drogue data points were available with the smaller L s /D 0 ratio: DDT-1, 
DDT-2, DDT-3, and a 10% scale Texas A&M University Wind Tunnel Test. 11 All of these tests were considered to 
be in a “clean” wake. The resulting distribution published in CPAS Model Memo vll did not have a central 
tendency, so a uniform dispersion was used. However, the additional data from HIVAS appears to show a central 
tendency, so a normal distribution was fit to the composite data, as shown in Figure 21. The EDU Pilot parachute 
also has half the number of gores as the CPAS Drogue, so some loss of drag efficiency was expected. 



30 35 40 45 50 55 60 

EDU Pilot Drag Area, C D S (ft 2 ) 

Figure 21. Comparison of expected Pilot parachute drag area with HIVAS axial flow results. 

VI. Operational Lessons Learned 

Some procedural issues were encountered during this test. The initial method to deploy the parachute into the 
airflow by pulling a line resulted in a premature deployment and HIVAS was shut down. A simpler method of 
manually tossing each canopy into the flow was successful. During the second run, the block holding the parachute 
began to move, despite being several times heavier than the force generated by the test parachute, and again the test 
was temporarily suspended. The parachute drag had caused a lever-arm effect on the block, tipping it enough to 
reduce the surface area until friction was no longer enough to hold it in place. This was mitigated by adding a 
second block to the test stand and securing the assembly to rails embedded to the deck with chains. 

VII. Conclusion 

The HIVAS facility at China Lake was used in lieu of flight test to determine parachute drag performance of two 
small CPAS parachutes. Another parachute with known performance was tested as a control. 

HIVAS produces a non-uniform flow, which was assumed to be axisymmetric for this analysis, though limited 
flow survey maps are available at certain settings and distances. Effective flow velocities at each throttle setting 
were computed for each size parachute using the definition of the momentum flux factor. The FBCP was not as 
stable as the other parachutes, requiring additional calculations of the projected flow at an off-axis angle. 

The computed drag coefficient on the control AFSAT Drogue was about 10% higher than previous data, which 
may indicate an under-estimation of the effective flow velocity. Steady-state data from HIVAS compared favorably 
to previously estimated parachute performance based on scaled data from analogous CPAS parachutes. Inclusion of 
the HIVAS data with previous data increases each standard deviation, making the model more likely to encompass 
actual performance in flight. The resulting drag area distributions will be included in the next version of the CPAS 
Model Memo. 
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The range of effective dynamic pressure tested was similar to the extent of the CPAS deployment envelopes, 
seen in Figure 22. However, HIVAS is limited by its ambient conditions in the Mach number and pressure altitude 
available. Larger parachutes cannot be tested effectively due to the flow decay from the center core and ground 
effects. Deployment was not representative to flight, so inflation parameters could not be obtained. 



Figure 22. Effective HIVAS flow velocity compared to CPAS deployment envelopes. 

Overall, HIVAS provided valuable data to the CPAS program with a much lower cost and complexity relative to 
flight testing. 
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